Open circular billiards and the Riemann hypothesis 
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A comparison of escape rates from one and from two holes in an experimental container (e.g. 
a laser trap) can be used to obtain information about the dynamics inside the container. If this 
dynamics is simple enough one can hope to obtain exact formulas. Here we obtain exact formulas 
for escape from a circular billiard with one and with two holes. The corresponding quantities 
are expressed as sums over zeroes of the Riemann zeta function. Thus we demonstrate a direct 
connection between recent experiments and a major unsolved problem in mathematics, the Riemann 
hypothesis. 
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Billiard systems, in which a point particle moves freely 
except for specular reflections from rigid walls, permit 
close connections between rigorous mathematics and ex- 
perimental physics. Very general physical situations, in 
which particles or waves are confined to cavities or other 
homogeneous regions, are related to well understood 
billiard dynamical systems, directly for particles and 
via semiclassical (short wavelength) theories for waves. 
Precise billiard experiments have used microwaves in 
metal 0, and superconducting |3l| cavities and with 
wave guides 0, visible light reflected from mirrors [3, 
phonons in quartz blocks electrons in semiconduc- 
tors confined by electric potentials H, |1| , and atoms in- 
teracting with laser beams |9| . llCt . Experiments em- 
ploy both two and three [13 dimensional geometries, 
and both closed and open systems. For instance, es- 
cape of cold atoms from a laser trap with a hole was 
studied in 0, 0|. Closed systems exhibit energy level 
distributions and scarring of wavefunctions as predicted 
by semiclassical |l3| and random matrix [l^ theories. 
Open systems exhibit unexpected and incompletely un- 
derstood phenomena such as fractal conductance fluctu- 
ations H, 0, 0, 01 . For both closed and open sys- 
tems, the behavior depends crucially on the classical dy- 
namics, which can be tuned to be integrable, chaotic, 
or mixed 0, 0, 0, 0|. Escape rate is a character- 
istic of open billiards which is both experimentally ac- 
cessible [3, and important for transport properties of 
many systems. In this paper we propose to ask what may 
be understood about the intrinsic dynamics of billiards 
using only this experimental escape information, specifi- 
cally by comparing systems with one and two holes. The 
two hole escape rate is not measured in the experimen- 
tal literature, but presents no fundamental difficulties. 
We naturally begin with the simplest of shapes, the cir- 
cle, and find remarkable exact expressions based on the 
most famous unsolved mathematical problem, the Rie- 
mann Hypothesis. 



At long times, the probability of a particle remaining 
in an integrable billiard with a hole is well known to ex- 
hibit power law decay, in contrast to exponential decay 
from strongly chaotic billiards 0|, however the coeffi- 
cient of the power ( "escape rate" ) in the integrable case 
has not been computed exactly to our knowledge. Nu- 
merical simulations can lead not only to uncertainty in 
the coefficient, but also misleading conclusions, for ex- 
ample the power law decay of an integrable billiard can 
be masked by an exponential term at short times. 

Here we consider the circle billiard, which is integrable 
due to angular momentum conservation. Some three di- 
mensional cases, namely the cylinder and sphere, can be 
treated analogously. The circle billiard of unit radius 
has collisions defined by the angle around the circum- 
ference (3 S (— TT, tt] and the angle between the outward 
trajectory and the normal -ip £ (— 7r/2, 7r/2). The billiard 
map is then (/?, ijj) (/5 + tt — 2^0, ip) where angles are 
taken modulo 27r as usual. The time between collisions 
is r = 2cosi/'- 

The dynamical evolution is of two types depending on 
the value of ip. For tp = ipm^n = 7r/2 — rmr/n with m < n 
coprime integers, the trajectory is periodic with period n 
and the (3 values are equally spaced at intervals of 27r/n. 
For Ip which are irrational multiples of tt, the trajectory is 
uniformly distributed in /? (see for example Prop. 4.2.1 
of Ref. ^]). 

For the escape problem the billiard is filled with a uni- 
form density of particles normalised to unity given by 
cos 4' dp dip / (Att) at the boundary; this is the usual equi- 
librium measure for two dimensional billiards. Note that 
for the circle, any function of ^ leads to an invariant mea- 
sure since ip is a constant of motion. Two (possibly over- 
lapping) holes are placed at the boundary at € [0, e] and 
/3 e [6,6 + e]; the one hole problem is simply 6 = 0. The 
number of collisions to escape is some function N(Pq, ipo) 
(possibly infinite) of the initial conditions and the time 
to escape is t = NT = 2cosipN{f]o,ipo). Considering 
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FIG. 1: Geometry of the billiard. 



the continuous and discrete time problems leads to es- 
cape times differing at most by the maximal free path 
of 2; only the discrete time calculation will be required 
for the leading (in time) behavior of escape, that is, the 
coefficient of 1/t for the density remaining at time t. The 
geometry of the billiard is shown in Fig. ^ 

From above, the only trajectories which do not escape 
are some of the periodic orbits with period n < 27r/e. 
For an arbitrary aperiodic (finite time) orbit we define 
the period n to be the smallest positive integer for which 
|0n — 0o| < £• For the largest possible value of n < 
27r/e, the n + 1 points (f>o, (/jit ■ ■ ■ 4'n automatically have 
a minimum distance < e on the circle, so the inequality 
is satisfied by this definition. In units of n collisions, the 
particle's orbit precesses slowly enough to be captured by 
one of the holes, so at very large values of t, only slightly 
precessing orbits remain. 

Now we compute the density remaining near periodic 
orbits (/3o, m, n) at long time t. Such a long lived trajec- 
tory has '0 = ^m.n + 'n foi' ?7 e. A prime will indicate 
values taken modulo 27r/n and lying in [0, 27r/n). Thus 
the dynamics is now /3' ^ /3' — 2r] and escape takes place 
when this passes one of the values e or e -I- 0' for 77 > 0, 
or or 0' for r] < Q. The set of initial conditions for 
which escape takes at least time t, hence t/(2 cos V'm.n) 
collisions is thus 
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for 77 > and similar expressions for 77 < 0. Integrat- 
ing over this region, including the cos "0 weighting for the 
invariant measure and a factor n to account for the equiv- 
alence modulo 27r/n, and summing over all non-escaping 
periodic orbits gives 
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where the sum is restricted to < m < n < 27r/e for 
coprime m, n, and faster decaying terms in t have been 
neglected (for example in approximating cos "0 by a con- 
stant near each periodic orbit). Note that n = 1 is in- 
cluded, but gives no contribution. Now we write the sin 
functions as sums of exponentials and apply the Ramanu- 
jan identity |2]| 
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where /i is the Mobius function, defined by /i(l) = 
1, /i(p) = —1 for primes p, ii{mn) = /i(TO)/i(n) if 
gcd(m,n) — 1, otherwise fj,(rnn) = 0. The result is 
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Here (j){n) is the Euler totient function, giving the number 
of positive integers m < n with gcd{m,n) ~ 1. Conven- 
tionally (/)(!) = 1. Of course Poo depends on both e and 
6, but this is suppressed for brevity. 

Poo is a finite sum, piecewise smooth when considered 
as a function of e and/or 9 with the number of terms 
unbounded in the limit of small e. For small holes we 
can extract the asymptotic behavior using Mellin trans- 
forms. Write P(s) — Pooe'*~^de and then Poo = 

2^ Ic~i^ ^~'^P{s)ds gives the required asymptotic se- 
ries as a sum of residues. The constant c is determined 
by convergence of the first integral. This integral is fi- 
nite (Poo = when e > tt) and so converges for suffi- 
ciently large s, leading to a constant c to the right of 
all poles of P(s). For sufficiently large s convergence 
of the sum is uniform in e, so interchanging the sum 
and the first integral, integrating over e and substitut- 
ing 9' = f {n9 / {2TT))2n / n where / indicates the fractional 
part, we find 
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Now we consider the case of rational angles 9 = 2'Kr/q 
where gcd{r,q) = 1; the r = 0, g = 1 case gives a single 
hole. In this case the sum over n splits into conjugacy 
classes modulo q and the fractional parts are known ratio- 
nal numbers which depend only on the conjugacy class. 
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Let us consider the required sum 

(p{n) - n{n) 
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by first dividing all quantities through by 6 = gcd(a,q), 
writing n' — n/h, a' — a/h and q' — q/h. We then have 



E 



n' =a' (mod q') 



bn') - ^l{bn') 



(7) 



with a' and q' coprime.Now we consider Dirichlet charac- 
ters [i^ x(?^')i defined as follows. The conjugacy classes 
modulo q' which are coprime to g' form an abelian group 
under multiplication, of order (jiW)- Since the group 
is abelian and finite, there are (t>{q') irreducible rep- 
resentations X which each n' coprime to q' is rep- 
resented by a complex root of unity xi^^') satisfying 
x(m')x(n') = xCt^'j^Oi xC^*') = if n' and q' have a 
common factor. 

A central result in representation theory is the orthog- 
onality theorem, here 



^^Xia')x{n') = Sa',n' 



(8) 



where 6a',n' = 1 if a' = n' (mod q'), zero otherwise, and 
the bar indicates complex conjugation. 

Inserting the orthogonality relation into the sum al- 
lows n' to be summed over all integers and decomposed 
into prime factors n' — YinP"''' 



E 



(j){n) - fi{n) 



n—a (mod q) 



T \i / n' = l ^ ' 

xK) ^n^^^^)""' <i>{^^)=m n (i-p"') 

^ [ m np(-l)"^ Ihn' for all c > 1 
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The condition on \x is taken into account by setting a^, = 
if p|6 and summing = 0, 1 otherwise. The sum for 
(/) is given by the a = term together with a geometric 
series. The result is 
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where characters are taken modulo g' and 



n— 1 p 



) (12) 



20 



10 



-10 



-20 



-20 



-10 



10 



20 



FIG. 2: Poles of P(s) (Eq.Ell for g = 1 (one hole) at -2, odd 
integers less than or equal to 1, and nontrivial values with real 
part —1/2 assuming the Riemann hypothesis. The contour of 
Eq.|K|is given. 



is the Dirichlet L function, which in the case q' = 1 (ie 
x{n) — 1 for all n) reduces to the Riemann zeta function. 
For each q' there is a trivial character for which x(a') ~ 1 
for all a' coprime to g'; this gives the same product as 
the Riemann zeta function, excluding primes dividing g', 
thus 
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Our first main result is the exact expression for the prob- 
ability P{t) of remaining in the unit circular billiard with 
two holes [0, e] and [27rr/g, iTxrjq -I- e] is thus 
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where, as above, b — gcd{a,q), a' — a/b, q' = q/b and 
the characters are taken modulo g'. In performing the 
contour integral (Fig. [SJ we recall that the contour lies 
to the right of all poles of the integrand (c > 1) and that a 
semicircular arc to the left may be added that avoids the 
poles and for which the integral vanishes in the limit of 
infinite radius. For large negative real part the integrand 
is controlled by the (27r)^"''^ and for large imaginary part 
the integrand is controlled by the cubic prefactor in the 
denominator. 

We now consider specific values of g; recall that the 
angle separating the holes is 6* = 27rr/g. Thus g = 1 is a 
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23(3+l){s+2){3» + l-l)C(3+l) 
2°(13C(3)+3°+^(C(s)-l)+3)-C(3)(3°+''+5) 
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[(7r/3) = + '-(6=+8.12''-25.30''+{l-3.2=-13.3=-8.4= 
+ 25.5= + 27.6°-25.10°+8.12=-25.15=+25.30'')C(s))] 
x[2s(s+l){s+2){2'' + l-l){3'' + l-l)f(s+l)]-l 



TABLE I; The function P(s)(g/27r)"+^ (see Eq.Eill for q 
1,2,3,4,6 and r = 1. 
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TABLE II: Some residues of P(s)e~'' given in table Q divided 
by the factor e"". The Ine appears for q = 6 due to a double 
pole at s = —1. There are also poles for further negative odd 
s, and along the critical line IZe s = —1/2. . 



single hole and q = 2 is two opposite holes. Experimen- 
talists could especially notice a contrast between small 
values of q in which the angles are simple rational multi- 
ples of TT and the escape rate is expressed in terms of only 
a few L-functions (in fact only the Riemann zeta func- 
tion for q = 1,2,3,4,6), and angles which are not close 
to rational numbers with small denominators. Particular 
values of P{s) are given in tableland residues of the triv- 
ial poles at integer s are given in table UTI Some points 
to note are that the leading behavior at s = 1, that is, 
of order is purely given by the total size of the holes. 
The q = 1 case is twice as large since there is only a 
single hole in this case. The next possible term would 
be s = 0, however the residue is zero since the C,{s + 1) 
in the denominator has a pole there. The s = — 1 term 
is more complicated for g = 6, involving a term in elne 
due to the double pole in this case. The s = —2 terms 
are remarkably simple, and are due to the s -I- 2 in the 
denominator. For these values of q the second to leading 
order terms, which are not given in the table, come from 
the nontrivial zeroes of the C,{s + 1), and for multiplic- 
ity m are of order •ye(lne)™~^ if all zeroes of C,{s) lie in 
TZe s < 1/2. The latter is a statement of the Riemann 
hypothesis jl^. An alternative formulation stating that 
all the nontrivial zeroes of C{s) have TZe s = 1/2 is easily 



shown to be equivalent using the functional equation |23| 
relating C(s) and C(l — s). Thus this celebrated unsolved 
problem is equivalent to either of 

lim lim e^-^^^{tPi{t) -2/e) = (15) 
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for every S > 0, and the subscript indicates the one or 
symmetric two hole problem. This is our second main 
result: The difference between the escape of the one and 
two hole problems is determined to leading order in the 
small hole limit by the Riemann Hypothesis. The gen- 
eralised Riemann hypothesis is the equivalent statement 
for L functions, and implies that corrections are of or- 
der e''+^/2 in the rational two hole case. For irrational 
9 the above analysis breaks down; presumably the poles 
on the critical line become dense, blocking analytic con- 
tinuation. However the leading order (in e) term can be 
shown to be 9 independent for 6* > 0, as for the ratio- 
nal values given above, as follows. For the leading order 
behavior, the sum over n can be approximated by an in- 
tegral, with parts of the summand replaced by "mean 
field" averages (). If 6* is irrational, the fractional parts 
are uniformly distributed, so that we compute 
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We also use {4){n)) = Qu/tt'^ and (/i(n)) = 0, so that 
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as required. 

We now present some numerical tests. The main re- 
sults of this paper are exact. However there are some 
related questions to do with rates of convergence of vari- 
ous limits which are of great importance to numerical or 
experimental extensions of this work. The rate of con- 
vergence of tP{t) as i ^ oo is considered in Fig. O The 
convergence of Ea. (|14|l with the number of residues is 
tested in Fig. 01 Note that considering only the real 
poles would lead to a Taylor series of integer powers, that 
is, an analytic function; the complex poles are essential 
to obtain the step structure. 

In conclusion, we computed the escape from a circular 
billiard with one hole and related it to the Riemann Hy- 
pothesis, computed escape for two holes (rational case) 
and related it to the generalized Riemann hypotheses and 
obtained the leading order behavior for the irrational two 
hole case. Many interesting questions remain, concerning 
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FIG. 3: The scaled survival probability etP{t) (see Eq. 
appears to approach a limiting function as e — > with et held 
constant. Here we use 10** random initial conditions, 6 = 1 
(an irrational multiple of tt) and the curves are e = 10~"^^ 
with n = . . . 6. At large et the function converges to unity, 
consistent with Eq. (I18II . 




FIG. 4: Numerical computation of P'^{e) for the single hole 
case q = 1, using Eq. (EU and varying the number of poles 
considered on the critical line; real poles are considered for 
s > —10 in all plots. From Q it is seen that the second 
derivative of Poo is a function with uniformly spaced steps 
and constant average gradient in the variables shown on the 
axes. Convergence is evident except for very small for 
which more real poles would be required; large require 
more critical poles for the steps to be visible. 



one or two hole escape from pseudo-integrable, chaotic or 
mixed bilhards, and to the quantum/wave signature of 
these systems 23] . A fuller understanding of open quan- 
tum billiards should have important practical benefits, 



for example in the design of microlasers |2J|. 

We would like to thank E. Croot and Z. Rudnick for 
helpful discussions. LB was partially supported by NSF 
grant DMS #0140165. CD is grateful for the hospitahty 
of the Center for Nonlinear Science, Georgia Tech. 



[4; 
[s: 
[6: 
[7: 

[9 
[10 

[11 

[12 



[is: 

[14 

[is: 



[16: 
[17: 



[18 
[19 

[20 

[21 

[22' 
[22, 

[24] 



S. Sridhar, Phys. Rev. Lett. 67, 785 (1991). 

J. Stein and H. J. Stockmann, Phys. Rev. Lett. 68, 2867 

(1992). 

H. D. Graf, H. L. Harney, H. Lengeler, C. H. Lewenkopf, 

C. Rangacharyulu, A. Richter, P. Schardt and H. A. Wei- 
denmuUer, Phys. Rev. Lett. 69, 1296 (1992). 

E. Persson, I. Rotter, H.-J. Stockmann, and M. Barth, 
Phys. Rev. Lett. 85, 2478 (2000). 

D. Sweet, B. W. Zeff, E. Ott, and D. P. Lathrop, Physica 
D 154, 207 (2001). 

C. EUegard, T. Guhr, K. Lindemann, J. Nygard, and M. 

Oxborrow, Phys. Rev. Lett. 77, 4918 (1996). 

T. Sakamoto, T. Takagaki, S. Takaoka, K. Gamo, K. 

Murase, and S. Namba, Japan J. Appl. Phys. 30, L1186 

(1991). 

J. P. Bird, J. Phys. Cond. Mat. 11, R413 (1999). 

V. Milner, J. L. Hanssen, W. C. Campbell and M. G. 

Raizen, Phys. Rev. Lett. 86, 1514 (2001). 

N. Friedman, A. Kaplan, D. Carasso, and N. Davidson, 

Phys. Rev. Lett. 86, 1518 (2001). 

M. F. Anderson, A. Kaplan, T. Griinzweig, and N. David- 
son, arXiv:quant-ph/0404118 (2004, unpublished). 
C. Dembowski, B. Dietz, H. D. Graf, A. Heine, T. Pa- 
penbrock, A. Richter, and C. Richter, Phys. Rev. Lett. 
89, 064101 (2002). 

J. Main, Phys. Rep. 316, 234 (1999). 

P. J. Forrester, N. C. Snaith, and J. J. M. Verbaarschot, 

J. Phys. A 36, Rl (2003). 

R. P. Taylor, R. Newbury, A. S. Sachrajda, Y. Feng, P. 
T. Coleridge, C. Dettmann, N. Zhu, H. Guo, A. Delage, 
P. J. Kelly, and Z. Wasilewski, Phys. Rev. Lett. 78 1952 
(1997). 

M. Terraneo, J. Mod. Opt. 49, 2027 (2002). 

R. Crook, C. G. Smith, A. C. Graham, I. Farrer, H. E. 

Beere, and D. A. Ritchie, Phys. Rev. Lett. 91, 246803 

(2003). 

L. A. Bunimovich, Chaos 11, 802 (2001). 

W. Bauer and G. F. Bertsch, Phys. Rev. Lett. 65, 2213 

(1990). 

A. Katok and B. Hasselblatt, Modem theory of dynamical 
systems (Cambridge University Press, 1995). 
H. Davenport, Multiplicative number theory (Springer- 
Verlag, New York, 1980). 

J. B. Conrey, Not. Amer. Math. Soc. 50, 341 (2003). 
R. W. Robinett and S. Heppelmann, Phys. Rev. A 65, 
062103 (2002). 

G. D. Chern, H. E. Tureci, A. D. Stone, R. K. Chang, 
M. Kneissel, N. M. Johnson, Appl. Phys. Lett. 83, 1710 
(2003). 



